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We consider the question of whether resonant inelastic x-ray scattering (RIXS) can be used to 
detect many-body excitations that are coupled to the spin-chirality terms Si • (Sj x S^) in a Mott 
insulator. We find that while the spin-chirality terms are in general absent in the usual experimental 
setups of RIXS, there are prospects of realizing such terms if one considers instead the scattering near 
a pre-edge. We then perform detailed analyses for the square and the kagome lattices, and brief 
analyses for the triangular and the honeycomb lattices, in which we show that the spin-chirality 
terms are indeed present in all the above lattices, but that they occur at a higher order in our 
expansion for the kagome and the honeycomb lattices. The merit of using RIXS in addition to 
Raman spectroscopy to detect excitations that are coupled to the spin-chirality terms is also briefly 
discussed in the context of the emergent gauge boson in the U{1) Dirac spin liquid. 
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I. INTRODUCTION 

With its ability to probe generic many-body excita- 
tions, Raman spectroscopy has become an important tool 
for understanding strongly correlated electronic systems, 
including the Mott insulators i Unfortunately, in Raman 
spectroscopy the photon momenta are generally negligi- 
ble when compared with the inverse lattice scale, thus 
making it essentially a zero-momentum-transfer probe 
and limiting its usefulness in certain cases. 

In the specific context of the U{1) Dirac spin liquid 
state in the spin-1/2 kagome lattice, in which the low- 
energy effective theory is described by chargeless spin-1/2 
fermions (spinous) coupled to an emergent C/(l) gauge 
field)^ we previously proposed Raman spectroscopy as 
a way to detect the spinon continuum and the fluctua- 
tions of the emergent gauge fieldi^ The prospect of de- 
tecting the emergent gauge field in experiments is partic- 
ularly significant, given the role such gauge fields have 
played in theories of quantum spin liquids. Unfortu- 
nately, our calculations also reveal that the signal coming 
from the emergent gauge field is suppressed by a factor 
of g^, where q is the momentum transferred to the sys- 
tem. Indeed, Raman spectroscopy on herbertsmithite 
ZnCu3(OH)6Cl2, a possible material realization of the 
[/(I) Dirac spin liquid state, has found a broad contin- 
uum in the spectrum that could be attributed to the 
spinon continuum, but has shown no signs of the emer- 
gent gauge field.— 

Given the large momentum carried by x-ray, one would 
imagine that resonant inelastic x-ray scattering (RIXS)^ 
may provide a better prospect of detecting the emergent 
gauge field. However, in our derivation, the detection of 
the emergent gauge field in Raman spectroscopy depends 
crucially on the coupling of the external photons to the 
spin-chirality terms • (Sj x Sfc) in the system, which in 
turn relies crucially on the link between the photon po- 
larizations and the direction of the virtual electron hops 



in the lattice induced by the virtual absorption and emis- 
sion of the photons. This link is in general absent in the 
current theoretical discussions and experimental setups 
of RIXS, in which the virtual absorption and emission of 
the photons are accompanied by intra-site electron hops. 



In this paper, we propose performing RIXS near a pre- 
edge, in which case the usual virtual processes with intra- 
site photon-induced electron hops are suppressed, thus 
allowing virtual processes with inter-site photon-induced 
electron hops to manifest. Indeed, inter-site dipolar con- 
tributions have previously been identified in the absorp- 
tion and Auger spectrum of^"— Ti02 and^iiii La2Cu04. 
To analyze the contributions by such inter-site processes 
to the RIXS signals, we modify the Shastry-Shraiman 
formalis m^ ^1^^ used in deriving the corresponding results 
in Raman spectroscopy, and show that the spin-chirality 
terms indeed appear in both the square lattice (cuprate) 
and the kagome lattice (herbertsmithite), but that the 
first appearance of such terms occurs at a higher order 
in the latter case. 



This paper is organized as follows: in Sec. |TT] the 
Shastry-Shraiman formalism in Raman spectroscopy is 
reviewed to set the stage for Sec. IIIIl in which the mod- 
ifications to this formalism to the case of RIXS are dis- 
cussed and illustrated. In Sec. IIVI the possibility of de- 
tecting the spin-chirality terms in RIXS is considered in 
this modified formalism and our new proposal is pre- 
sented, followed by detailed analyses for the square and 
the kagome lattices, as well as brief discussions for the 
triangular and the honeycomb lattices. Further discus- 
sions ensue in Sec. |Vl in which our motivating example 
of the U{1) Dirac spin liquid is considered again to put 
our proposal into perspective. 
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II. REVIEW OF THE SHASTRY-SHRAIMAN 
FORMALISM IN RAMAN SPECTROSCOPY 

In the Shastry-Shraiman formalism, the electron- 
photon interaction He is treated as a time-dependent 
perturbation on the time- independent Hamiltonian i?ind- 
The latter consists of the Hubbard Hamiltonian i/nb = 



T-- U 



tonian H 



ni-^riii^ and the free photon Hamil- 
q (here a labels the photon 



7 — Z^q,a'^q"'q 

polarizations). Applying Fermi's golden rule, the transi- 
tion rate Wfi from an initial state |i) to a final state |/) 
is given by: 



Wf, = 2i,\{J\T\i)\H{£f^£^^ 



(1) 



where T is the scattering T-matrix. Keeping only terms 
that are second order in the photon operators , the T- 
matrix can be decomposed as T = Tr + Tnr, with Tr the 
resonant part and Tnr the non-resonant part, of which 
only the former is important in a Mott insulator. Next, 
the energy denominator that appears in Tr is further ex- 
panded by treating the hopping part Ht = J2ij,a ^ij^la^j-y 
of the Hubbard Hamiltonian as a perturbation on the re- 
maining terms in Hind- To be more precise. 
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FIG. 1: Virtual processes that contribute to the resonant Ra- 
man scattering T-matrix Tr. The process depicted in (a) con- 
tributes to the Fleury-Loudon term while the ones depicted 
in (b) and (c) contribute (among others) to the spin-chirahty 
terms at the leading order. Here and henceforth thiclc (blue) 
arrows denote electron hops that are accompanied by virtual 
absorptions or emissions of photons, thin (magenta) unbro- 
ken arrows denote movements of electrons in non-photon- 
induced internal hops, and thin (magenta) broken arrows de- 
note movements of holes in non-photon-induced internal hops. 
The order of hops is indicated by lowercase roman letters next 
to the corresponding arrows. 



Tr ^H, 



(1) 



1 



£i - [Huh + Hj) + i-q 



-H. 



(1) 



£i — Hjj — Hj + irj 
^[ '£,-Hu-H^ + iV 



H, 



c ' 



(2) 



(3) 



where £i is the initial energy of the unperturbed sys- 
tem, H^^ consists of the terms in He that are first order 
in the photon operators Cq, and Hu = C/X)i"it"i-1 
the on-site Coulomb repulsion in the Hubbard Hamil- 
tonian. In the context of Raman scattering, H^^ ^ 
'Etj,a^tvcLcj<yigiec,a'^e'^'' + gfepa^e-''^''), where gi 
(gf) denotes the appropriate coupling constants between 
the electron and the incoming (outgoing) photon, x = 
(xi -|-Xj)/2 is the mid-point between site i and j, and 
(e^) and k (q) denote the polarization vector and mo- 
mentum of the incoming (outgoing) photon, respectively. 

In a Mott insulator, Tr connects between two spin 
states (i.e., states with zero double occupancy) and hence 
can in principle be expressed in terms of spin operators. 
In the Shastry-Shraiman formalism, this is achieved by 
inserting a complete set of states (in the lattice occu- 
pation basis with a fixed spin quantization axis) in be- 
tween the operators in Eq. ([3]), under which the energy 
denominators {£i — Hu — H-y+iri)~^ become c-numbers.— 
Moreover, once an initial spin state (in that same basis) 
is specified and a particular choice of individual term is 
chosen for each H^ and h''^^ in Eq. ([3]), the intermediate 
states are uniquely determined and thus can be trivially 
rcsummed. Hence the matrix elements of Tr with respect 



to spin states can be expressed as a sum of chains of elec- 
tronic operators. These chains of electronic operators can 
be visualized as virtual processes in which electrons hop 
around the lattice and can be converted into spin oper- 
ators by the identities Xaa' = c^'Ccr = + S • t„„i 



For t <ti U and near 



and Xcrcr' = c<tC^ 
resonance (i.e., uji w U , in which cji is the energy of the 
incoming photon), the contributions to Tr are dominated 
by virtual processes in which all intermediate states have 
exactly one hole and one doubly occupied site (a.k.a. dou- 
blon). In such case {£i — Hu — Hj + irj) sa (wj — U), and 
the matrix elements of Tr can thus be organized as an 
expansion in t/^Ui — U). 

To the lowest nontrivial order in t/(u)i — U), the 
T-matrix obtained in the Shastry-Shraiman formal- 
ism reproduces the Fleury-Loudon Hamiltonian iJpL = 

Er.r' u^i^i • • m)(1/4 - Sr • Sr')^ and is con- 

tributed by virtual processes of the form shown in 
Fig. 1(a) At the f^/ {uji — U)^ order, individual processes 
that contribute to the spin-chirality terms, such as the 
ones shown in Figs. 1(b) and 1(c) start to appear. How- 



ever, in the case where only nearest-neighbor hoppings 
are included, the sum of their contributions is found to 
vanish in the square and the triangular lattices, while 
it remains nonzero in the honeycomb and the kagome 
lattices^ It is worth noting that the spin-chirality terms 



appear exclusively in the (e^e^- 
ncl in the scattering T-matrix. 



e^ef ) polarization chan- 
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III. MODIFICATION OF THE 
SHASTRY-SHRAIMAN FORMALISM TO RIXS 

Since both Raman scattering and RIXS are resonant 
two-photon processes, it should be possible to modify the 
Shastry-Shraiman formalism to the case of RIXS. Indeed, 
similar expansion of Tr had been made by van den Brink 
and van Veenendaal.— However, in that work the denom- 
inator in Eq. ^ was expanded by treating the (appro- 
priately modified, see below) Hubbard Hamiltonian as 
a perturbation on the free photon and the atomic (see 
below) Hamiltonians. In practice, the Hubbard Hamil- 
tonian is then expanded in the usual way one derives 
the Heisenberg Hamiltonian i ^^'^^ For the present work, 
we follow instead the spirit of the Shastry-Shraiman for- 
malism and take the terms in i?ind in which site indices 
change as perturbations on the free photon and the on- 
site terms in i?ind- Since the main purpose of this paper 
is to identify virtual processes that may give rise to the 
spin-chirality terms, of which the two expansion schemes 
agree except that the prefactors coming from the energy 
denominators are organized differently, we shall not dwell 
on the relative merit of these two expansions, which may 
depend on one's identification of the resonant energy. In- 
stead, we simply state here the necessary modifications 
to the Shastry-Shraiman formalism in the case of RIXS. 

First, -ffp ^ now corresponds to virtual transitions in 
which an electron hops from a core state to a valence 
state while a photon is absorbed, or in which an elec- 
tron hops from a valence state to the core state while 
a photon is emitted. Hence we may write iJ^"* — 
'Ecv igzeaclJ^^Cca'^e'^''+igfepcl{j'')^^Cva^e''^''\wheTe 
the subscript c {v) labels a core (valence) state and J is 
a (possibly polarization-dependent) matrix that accounts 
for the matrix elements of the atomic transitions^^ As 
usual it is necessary to include only the core and valence 
states that are near resonance in 

Second, iJind raust now include extra terms that ac- 
count for the single-particle energies of the core states 
and possibly of the high-energy valence states, as well 
as their interactions with the low-energy valence states. 
Schematically, we can write i/ind = +i/atomic +^Hbi 
in which the atomic Hamiltonian iJatomic accounts for the 
energy difference between the core state and the valence 
state excited from it, while the modified Hubbard Hamil- 
tonian accounts for the interactions of the valence 
electrons among themselves, with the lattice potential, 
and with the core hole. In particular, since Hu captures 
only the low-energy effective Coulomb repulsion among 
the low-energy valence electrons, in principle it is nec- 
essary to include in generic Coulomb interactions 
c'pC^cs in which at least one of the electron oper- 
ators corresponds to a core or high-energy valence state. 
However, on physical ground it may be argued that the 
dominant effect of the core or high-energy valence state 
on the low-energy valence states would be modifications 
to the hopping parameters and the on-site potentials, i.e.. 
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FIG. 2: The lowest-order virtual processes that contribute to 
the one-magnon excitation in 2p — ^ 3d RIXS. 



iccVc,ijclcj and clceVe,ijclcj, where 



the subscript c (e) labels a core (high-energy valence) 
state while i,j are site labels of low-energy valence states. 
In practice, the effect of the core hole may be well cap- 
tured by an on-site potential Uc localized at the site 
where the core hole is present 

Third, unlike in Raman scattering, the core hole in 
RIXS has a very short lifetime and can decay via Auger 
processes. This introduces an uncertainty to the core- 
hole energy, which can be captured by replacing the in- 
finitesimal rj in Eqs. ([2]) and ([3]) by a finite energy broad- 
ening r. 

To illustrate the modified Shastry-Shraiman formal- 
ism, we briefly outline the virtual processes that con- 
tribute to the single-magnon excitation in 2p — > 3d 
RIXS and to the two-magnon excitation in Is — )■ Ap 
RIXS in cuprates, both of which have previously been 
proposedi^^ and experimentally studied.— 

For the 2p — )■ 3d RIXS, the lowest-order processes that 
contribute to the single-magnon excitations are purely 
atomic in nature and involve simply the photon-induced 
virtual transitions of a 2p core electron to and from the 
3d valence states (see Fig. [2] for illustration). Since the 2p 
states are spin-orbit coupled, spin is not a good quantum 
number, and a spin flip of the 3d electrons can occur, as 
long as it is accompanied by an appropriate change in the 
photon polarization. In the modifled Shastry-Shraiman 
formalism, the chain of electron operators that is associ- 
ated with the process depicted in Fig. [2] is: 



oc 



^ ] ('^p('^^)p(T^o-)(Cct'<^(t'p'Cp') 



tr{Af"'^x} 
m"^ - 2m"'3 • S 



(4) 



where the subscript p labels the 2p states while the 
subscript a labels the spin of the valence 3d states. 
Also, Xcrcr' = Co-Cct' ~ ^^cru' — S • Tg-cr' as bcforc while 
M:P^ = mfS.,^ + . r.,. = EpJ^iJ^C- Note 
that we have adopted a matrix convention for spin in- 
dices in the second line (which will be similarly adopted 
henceforth). From Eq. it can be seen that the spin- 
flip term arises from structure of the 2p — > 3d transition 
matrix elements. 
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■ 4p 

• 3d 



■ ip 



-S- 3d 



(a) 



(b) 



FIG. 3: The lowest-order virtual processes that contribute to 
the two-magnon excitation in Is — s- 4p RIXS. 



For the Is — >■ 4p RIXS, two-magnon excitation occurs 
when the core hole and 4p valence electron "shake up" the 
low-energy valence electrons and induce a pair exchange. 
Two such processes are depicted in Fig. [3l The chains of 
electron operators that are associated with the processes 
depicted in Figs. |3(a)] and [3(b)] are, respectively: 

^2-magnon i'^ls i"^^ )sp'^ip)i^ij'4'^j ) i^ji'^l'^i)(.'^ip' Jp's'^^is' ) 

= tr{J"(jt)^}tyt,, tr{xJX^} 



(5) 



^2-magnon ^ (^L ("^^ )fp'^ip) (''jicici) (i^j cj^Cj ) (c|^p, Jp, ^.z Ci^' ) 

= tr{J"(jt)^}t,,t,,tr{x.X,} 



(6) 



where N"^ = tr{ J"( jt)^}, the subscript s (p) labels the 
Is (4p) states, and the electron operators with no orbital 
labels are assumed to be that of the valence 3d states. For 
brevity here and henceforth we omit the spin indices in 
virtual hops that involve only the electrons, assuming 
that they arc appropriately summed within parentheses. 



Thus, 



c]„Ci 



Similarly, here and 



henceforth the sums over the spin and (for core hole and 
high-energy valence electrons) the orbital indices arc as- 
sumed in the two photon- induced hops, such that, e.g., 
{4p'Jp's'C^s') = T,p',s'4p'Jp's'C^s'■ Note also that the 
matrix notation in the second line of Eqs. ([5]) and ([6]) is 
now extended to include the orbital indices of the core 
and the 4p electrons. 

Under the assumption that the only effect of the Is 
core hole and the 4p valence electron is to introduce 
an additional potential Uc at the site i of which the 
atomic transition occurs, the coefficients of ri and 

^2-'magnon ^^lat comc from the energy denominators are, 
respectively: 



a 



a 



2-mag„on [^Suj + iT J SiJ - {U + Uc) + iV ' 

(8) 



(6) 



2-mag„o„ \^SuJ + iTj dcO - (U - Uc) + iT 




(a) 



(b) 



FIG. 4: Two processes whose contributions to the spin- 
chirality terms cancel out each other. 



Here Slj = LUi — {£4p—£is) is the detuning from the atomic 
transition. We remark that one must subtract from 
C2-magnon the Corresponding coefficients with Uc = to 
obtain the actual contributions of these two processes 
to the two-magnon transition, since the contributions 
of these two processes with Uc set to 0, together with 
other virtual processes in which the intermediate virtual 
exchange docs not involve the site i, constitute a part 
of the T-matrix that is proportional to the Heisenberg 
Hamiltonian and thus is not responsible for any actual 
transitions. 



IV. SPIN-CHIRALITY TERMS IN RIXS 

To obtain the spin-chirality terms, at least three lat- 
tice sites must be involved in the virtual processes. For 
example, one might consider the higher order processes 
in the Is — 4p RIXS shown in Figs. 4(a) and |4(b)] From 
the modified Shastry-Shraiman formalism, the associ- 
ated chains of electron operators are, respectively: 



TtLs^i4siJ^)lp^^p)iUe4c,){te,clc,^^^^^ 

= tr{J"{J^)^}tutiJtJ,tr{xiXJX^} 

= N^Hatejt,, (2iSf • (Sj X SO + . . .) , (9) 

T^%.s^i4siJ%c.p)ikj4cj)ik^4c,){^^^^^^ 

= tT{r{J^)^}t,jt,itu tT{xJXeX^} 
= N°'hutijtj^ {2iSj ■ {Si X S, 



.) + •■■) 



(10) 



Moreover, it can be seen that the coefficient coming from 
the energy denominators is Cs-sites = (<5w -I- iT)^'^(5u) — 
[U — Uc) + ir)~^ for both processes. Hence, while each 
process by itself contributes to the spin-chirality terms, 
the sum of their contributions vanishes Similar cal- 
culations show that the cancellation also occurs in the 
two analogous virtual processes in which electrons hop 
around the four-site loop in the square lattice. Further- 
more, it can be shown that such cancellations also occur 
for similar virtual processes in the 2p — )• id RIXS. 
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In Raman scattering, the analogous processes, in which 
the hops i and t; in Fig.[4]do not exist and which the hops 
ii and iv are photon induced, do not cancel out each 
other. Instead, the anticlockwise loop in Fig. 4(a) 



con- 



tributes to the e^ef photon polarization channel while 
the clockwise loop in Fig. |4(b)| contributes to the ejef 
channel, resulting in a nonvanishing contribution to the 



(e 



e^ef ) channel when summed. This suggests that 



in order for the spin-chirality terms to be realized in the 
scattering T-matrix, it is crucial for the photon polar- 
izations to be coupled with the directions of inter-site 
electron hops — a link that does not appear in the usual 
RIXS setups. 

That said, one should also note that the dipole mo- 
ment between a core orbital at one site and a valence 
orbital at one of its nearest-neighbor site is in general 
nonvanishing. Thus, in principle, a photon from the inci- 
dent x-ray beam can also induce a core-to- valence excita- 
tion across the two sites. Such an inter-site transition is 
in general suppressed by the reduced wavefunction over- 
lap and is thus masked by the corresponding intra-site 
transition. Moreover, for hard x-ray the distance be- 
tween two nearest-neighbor sites may also be equal to 
many wavelengths of the incident x-ray, which further 
reduces the transition amplitude for such an inter-site 
transition at near-horizontal incidence (relative to the 
two-dimensional lattice plane) as a result of the rapid 
oscillation of the electric field across the two sites. 

However, if the frequency of the incident x-ray is tuned 
to that of a forbidden atomic transition (e.g., the Is — >■ 3d 
transition in the Cu^"*" materials), then the near- resonant 
dipole inter-site transition needs only to compete with a 
near-rcsonant quadruple intra-site transition and a de- 
tuned dipole intra-site transition. With adequate lu- 
minosity, this may allow the signals from the inter-site 
transition to manifest in the spectrum. Indeed, contri- 
butions from such inter-site transitions have previously 
been identified in the x-ray absorption and Auger spec- 
troscopy oi^-^ TiOs and^iia La2Cu04. 

Moreover, for nearly two-dimensional materials such as 
cuprates and herbertsmithite, the rapid oscillation of the 
electric field between two lattice sites can be alleviated 
by arranging the x-ray to be at near-normal incidence 
relative to the two-dimensional lattice plane. Generally, 
by tuning the angle of incidence, the electric field across 
two nearest-neighbor sites can be made relatively uni- 
form while a sufficiently large in-plane momentum of the 
photon is maintained, such that a significant portion of 
the Brillouin zone can be explored. Furthermore, one 
may also consider resonances induced by soft x-ray (e.g., 
the 2s — )■ 3o? and the 3s — >■ 3d resonances in cuprates), 
which have larger wavefunction overlaps between the core 
orbitals and their nearest-neighbor valence orbitals. 

For the rest of this section we shall assume that the 
signals from such inter-site transitions can indeed be de- 
tected and consider in detail whether the spin-chirality 
terms can indeed arise from such a case. Specifically, 
we shall focus on s — 3c? inter-site transitions in Cu^+ 





FIG. 5: Orientations of the 3d orbitals in (a) cuprates and (b) 
herbertsmithite, and their effects on the signs of the photon- 
induced hopping amplitudes in the s 3d RIXS. Here the red 
(solid) and cyan (shaded) fillings indicate the relative signs of 
the angular part of the electron wavefunctions. Note from the 
figure that none of the photon-induced hopping magnitude is 
required to vanish by symmetry. 



materials with the square and the kagome lattice geome- 
tries, having in mind the realistic materials of cuprates 
and herbertsmithite. We shall also briefly comment on 
the cases of the triangular and the honeycomb lattices, 
in which the derivations of the spin-chirality terms are 
closely related to that of the square and the kagome lat- 
tices, respectively, and in which the former may be rel- 
evant to the new spin liquid candidate^^ Ba3CuSb2 09. 
For brevity we shall drop the factors gi and gf that are 
common to all virtual processes. 

In such virtual processes with photon-induced inter- 
site hopping, it is easy to check that the intermediate 
state obtained after a photon-induced hop has an en- 
ergy denominator of £]j = duj — (U — Uc) -\- iT, where 
Su! = u!i — {E^d — ^^corc) IS again the detuning from the 
atomic transition. The resonant condition is thus given 
by uji « {Sad- Score) + (U -Uc), under which « iT. For 
cuprates, U « 8.8 eV and Uc « 7.0 eV, while F « 0.75 eV 
for the K-edgei^ In comparison, in the effective one-band 
Hubbard model for cuprates, t sa 0.4 eV.— Observe that 
the expression of £d involves Uc, suggesting that the true 
resonant frequency of the inter-site dipole transition is 
offset from that of the intra-site quadruple transition by 
Uc- Indeed, the relative frequency shift between the inter- 
site dipole transition and the intra-site quadruple transi- 
tion has been used to explain the "three peaks" feature 
of Ti pre- if -edge absorption spectra in Ti02i^ The exis- 
tence of such frequency shift may thus allow the signals 
from the intra-site quadruple transition to be further sup- 
pressed relative to inter-site dipole transition when the 
frequency of the incident photon is tuned. 

If wc further assume that T U and that U and Uc 
arc of the same order of magnitude, assumptions that 
appear to be valid in cuprates, then the virtual processes 
are dominated by those having intermediate states with 
exactly one low-energy valence doublon not located at the 
core-hole site, and we can organize the T-matrix as an 
expansion in t/T similar to that in the Raman case. Since 
the effect of the photon polarizations is now mostly re- 
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tions that contribute to the spin-chirality terms in the s ^ 3d 

RIXS. FIG. 7: Mapping of two-internal-hop processes between the 

square and the triangular lattices. 



fleeted in the directions of the induced electron hops and 
since the spin-orbit coupling is negligible, we may take 



H, 



for the electron hop associated with the virtual absorp- 



tion of photon and ' ^ ± J(eQ • jjb) 

±J(eQ • m)(c|-i-^.s'^s) foi' tl^s electron hop associated with 
the virtual emission of photon, in which operators with 
the subscript s correspond to the s orbitals while opera- 
tors without orbital labels correspond to the valence Zd 
orbitals. Note that J is now a real scalar constant. The 
± signs in the above equations are determined by the rel- 
ative orientations of the d orbitals and are illustrated^^ 
in Fig. [SJ 

Now consider the particular case of the square lattice 
with only uniform nearest-neighbor hopping t. For such 
a lattice, virtual processes that involve valence electrons 
on at least three sites first appear at the order of two 
internal hops (i.e., at the (t^ J^ /£'n) -th order). Two such 
processes arc depicted in Figs. 6(a) and |6(b)| The corre- 
sponding contributions to the T-matrix are, respectively: 



rp{a) 

sq 



7^(e/^cl ,,C3)(t4c4)(i4c2)(e^j4ci,s) 
-j^-e^e^ tr{x3X4X2} 

'-'D 



S3 • (S4 X S2) 



■-D 



(11) 



= ^(-e^Jctc2)(t4c4)(teic3)(-e:f j4ci,.) 



■e/ef tr{x2X4X3} 



-^t^J^ie^ef) S2 • (S4 X S3) 



-D 



(12) 



tion of the spin-chirality terms have been adopted on the 
fourth line of Eqs. (HI]) and (fT2|l . 

For a fixed core-hole site, at this order there are three 
additional pairs of processes that contribute to the spin- 
chirality terms, which can be obtained from the processes 
depicted in Fig. [6] by successive 90° rotations about the 
core-hole site. Summing over all these processes, to this 
order the contribution to the spin-chirality terms by a 
core hole at site i is given by: 



rj-il 



(^i 



+ 



(13) 



Summing over all possible core-hole sites, and now restor- 
ing the exponential factor e**^'''"*'^ we havei^ 



E 

R 
X 



2ie.p 



+ 



R 



R 



(14) 



where the subscript R next to the parenthesis labels the 
site with which the spin-chirality terms are associated. 
From Eq. (fT4|) we see that for the square lattice there are 
indeed contributions to the T-matrix that couple to the 
spin-chirality terms at momentum q = — (i.e., the 
momentum transferred by the photon). 

The above analysis for the square lattice can be readily 
modified to the case of the triangular lattice, since pro- 
cesses with less than two internal hops can involve valence 
electrons at at most two sites and hence do not give rise 
to any spin-chirality terms, while the two-internal-hop 
processes in the triangular lattice and the square lattice 
are topologically the same (see Fig. [7] for illustration). 
For example, the contribution to the spin-chirality terms 
by the process depicted in Fig. [7Jb) can be read off as: 



(b) 



2it'J.Jf 



£■3 



(15) 



where = denotes the part of the T-matrix that contains 
the spin-chirality terms, and that a graphic representa- 



where the superscript w in the photon polarization ej 
corresponds to the unit vector w as depicted in the figure. 
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FIG. 8: Topologically distinct resonant virtual processes with 
inter-site photon-induced transitions in the s — ^ Sd RIXS in 
the kagome lattice, up to three internal hops. 



Assuming that all photon-induced hops have the same 
amplitude, so that Ji = Jf = J for all processes, we can 
sum up all contributions as in the square lattice case to 
obtain: 



2£l 




Hence, as in the square lattice case, the spin-chirality 
terms do appear in the triangular lattice at the 
[t^ .P I ED-th order. 

Next wc consider the kagome lattice with only uniform 
nearest-neighbor hopping t. In Fig. [8]wc list all the dis- 
tinct topo/o(7ies (as opposed to geometries, such that, e.g.. 
Fig. 8(c) is also representative of processes in which the 
third site is located at other locations) of resonant vir- 
tual processes up to three internal hops. Unfortunately, 
none of these processes generate any spin-chirality terms. 



FIG. 9: Virtual processes with inter-site photon-induced 
transitions that contribute to the spin-chirality terms in the 
s — >■ 3d RIXS in the kagome lattice. 



To see this, first observe that the processes depicted in 
(a) I |8(d)| involve valence electrons at only two 



Figs. 

fewer sites and thus cannot generate any spin-chirality 
terms (note that only core electrons are involved at site 
1 in all of the processes listed in Fig. [8]). Next, to rule out 
the processes depicted in Figs. 8(e) 8(g) note that a spin 



state is annihilated by two successive creation operators 
or two successive annihilation operators on the same site, 
regardless of the spin characters of these two operators. 
Consequently, if a site is transversed more than once in a 
virtual process, then each internal loop contributes to a 
separate trace in the Shastry-Shraiman derivation. For 
instance, corresponding to Fig. 8(g) we have: 



T^ifLe = e^/e- S^(cLc2)(4c30(4'Ci')(4'C2)(4ci,) 



3 72 



-X X ^ J 

£ 



■(cl.,,C2)tr{x3'Xi'}(4ci,.'5) 



t^J^ 

6/ef^^tr{X2}tr{x3'Xi'} 



(17) 



Similarly, T^^'^^^ 



(1/2- 



282-83) and so doesTHl^^. 



For the kagome lattice, the spin-chirality terms first ap- 
pear at the (i^ J^/£|))-th order, which arise from virtual 
processes in which the doublon hop through a hexagon. 
One such process is depicted in Fig. 9(a) whose contri- 
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bution is given by: 



T^^L = e^ef ^(cl,,C6)(4c5)(4c4)(cic3)(4c2)(4ci.. 

= 6/ tr{x6X5X4X3X2} 



^f^i Of 5 



D 

2£ 



J2 S„ • (Sb X SJ , (18) 

^ 6>a>b>c>2 

where the superscript w in corresponds to the unit 
vector w as depicted in the figure. One can check that 
processes with the same topology as the one depicted 
in Fig. 9(a) sum to a nonzero contribution to the spin- 
chirality terms in the {e^ef — e^ef ) channel at momen- 
tum q equal to that transferred lay the photon, which for 
brevity we shall not write down explicitly. It can also be 
checked that any process at this order with a different 
topology does not contribute to any spin-chirality terms. 

Compared with the corresponding terms in the square 
lattice, the spin-chirality terms in the kagome lattice are 
down by a factor of (t/fu)^, which can be significant in 
the limit where t <^ T even if the resonant condition is 
met. In such case one may want to consider also pro- 
cesses in which not all energy denominators are equal 
to £d- With this relaxed criterion, contributions to the 
spin-chirality terms can be found at the order of two in- 
ternal hops, in which the doublon hops through the core- 
hole site (see Figs. [9(b)] and [9(c) for illustrations). For 
instance, the process depicted in Fig. [9(b)] contributes: 

^2 j2 

/ < 72TF— ^TTT (^L C3 ) ( 4 ci ) (4 C2 ) (4 ci , , ) 



(C/c 



kagome 



tr{X3XiX2} 



£1{£d + Uc) 



S3 • (Si X S2) 



(19) 



Again it can be checked that all processes with the same 
topology as the one depicted in Fi g. [9(b)[ (which in- 
cludes the one depicted in Fig. 9(c) I sum to a nonzero 
contribution to the spin-chirality terms in the {e^ef — 
e^ef) channel. It is worth noting that such "back- 
tracking" processes are also present in the square lat- 
tice and carry opposite signs from the ordinary ones 
depicted in Fig. [6l Thus in the limit where t <^ T, 
the ratio of prefactors in the spin-chirality terms in the 
kagome lattice over that in the square lattice is given by 
{£d + C/c)"V i^B' - i^D + U,)-') = £d/U,. It is, how- 
ever, worth noting that t/T is not expected to be small 
in the case of cuprates. 

Since the honeycomb lattice has the same hexagon 
loops as in the kagome lattice and has no shorter (in 
terms of the number of hops) loops, it can be readily 
checked that the spin-chirality terms again first appear 
in the honeycomb lattice at the (i^ J^/£|))-th order when 
t <\£d\ and at the {t^P/{£D + uj)£n)-t\y order when 



t <^ \£d\, with Figs. 9(a) and 9(c) the typical contribut- 
ing processes in the respective cases. 



V. DISCUSSIONS AND CONCLUSIONS 

In this paper, we consider the question of whether 
RIXS can be used to detect many-body excitations that 
are coupled to the spin-chirality terms in a Mott insula- 
tor. We find that the spin-chirality terms are in general 
absent in the usual experimental setups, in which the 
spectroscopy is done near an absorption edge. The ab- 
sence of the spin-chirality terms in these setups can be 
traced to the lack of linkage between the virtual elec- 
tron hops and the photon polarizations. However, we 
argue that RIXS still holds a prospect of observing the 
effects of the spin-chirality terms if one instead considers 
spectroscopy near a pre-edge, in which case the intra-site 
dipole transitions are forbidden. 

Focusing on the Cu^"*" materials with the square and 
the kagome lattice geometries, we find that the spin- 
chirality terms are indeed presented in both cases under 
our new proposal. However, in the kagome case such 
terms appear only at a higher order in our expansion. 
In addition, we also find that as far as the spin-chirality 
terms are concerned, the scenario for the triangular lat- 
tice is analogous to that of the square lattice, while the 
scenario for the honeycomb lattice is analogous to that of 
the kagome lattice. It is worth noting that the situation 
we encounter in RIXS is essentially the reverse of what 
happens in the Raman case, in which the spin-chirality 
terms occur at the (t^/J7^)-th order in the kagome and 
the honeycomb lattices but not in the square or the tri- 
angular lattices. 

In comparison to the similar scheme to detect the spin- 
chirality terms in Raman spectroscopy ; ^ ^1 which had al- 
ready been realized)^ the present scheme in RIXS suffers 
from the reduced wavefunction overlaps in the inter-site 
dipole transitions. However, it has the advantage that 
excitations with finite momentum can be probed. To 
put this into perspective, let us return to the motivation 
we presented in the introduction, namely the emergent 
gauge boson in the U(l) Dirac spin liquid. In the U{1) 
Dirac spin liquid, the spin-chirality terms in the T-matrix 
correspond to flux- flux correlators, viz.: 



^M^/.=^2^|(/|T|^)p<5(^:/-f,0 

/ / 

- (i|6(Ak,Aw)6(O,0)|i) + •• 
q2e(Aw - VF^k) 

(Acj2 - „2 Afc2)l/2 + 



(20) 



(21) 



where Q denotes the step function, vp is the Fermi ve- 
locity at the Dirac cone of the U{1) Dirac spin liquid, 
b is the "magnetic field" associated with the emergent 
gauge boson, and Aw = — w/ (Ak = — k/) is the 
energy (momentum) transferred from the photon. If we 
assume that £0 in RIXS and (w^ — U) in Raman spec- 
troscopy are of the same order, the intensity of the sig- 
nal from the gauge boson in RIXS will be modified from 
that in Raman spectroscopy by a factor roughly equal to 
J^£d li^(£D + C/c) or J"^ j£\i, depending on which limit 
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one considers in RIXS. However, such comparison is not 
particularly meaningful since we have not considered how 
the background signals compare in the two cases. How- 
ever, the advantage offered by RIXS is not so much in the 
intensity of the signal but rather in its lineshape. In the 
Raman case where Ak « 0, the signature of the emergent 
gauge boson can manifest only as a power-law behavior 
near zero energy transfer, which can easily be masked 
by the elastic or quasielastic peak. In contrast, in RIXS 
the signal from the gauge boson has a sharp threshold at 
Auj = fpAfc, which varies as Ak varies. Thus, assuming 
modest intensities of the signals, it would be much easier 
to discern the emergent gauge boson in the case of RIXS. 

Of course, one should not underestimate the experi- 
mental challenges in realizing the proposal laid out in 
this paper. However, enormous progress in RIXS has 



been made in recent decadeS)^ with two-magnon excita- 
tions being observed^^^ and three-magnon excitations 
being proposed.— It is our hope that our proposal will 
further stimulate new theoretical and experimental ad- 
vances in the field. 
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